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LONG  WAVELENGTH  LIMIT  OF  THE 
CURRENT  CONVECTIVE  INSTABILITY 


I.  INTRODUCTION 

The  current  convective  instability  has  recently  been  suggested  as  a 
mechanism  to  generate  density  irregularities  in  the  auroral  ionosphere 
(Ossakow  and  Chaturvedi,  1979;  Vickrey  et  al.  ,  1980;  Chaturvedi  and  Ossakow, 
1981;  Keskinen  and  Ossakow,  1982;  Xeskinen  and  Ossakow,  1983).  These 
irregularities  can  cause  the  scintillation  phenomena  observed  by  the  DNA 
Wideband  satellite  during  periods  of  diffuse  aurora  (Fremouw  et  al.,  1977; 
Rino  et  al.  ,  1978),  and  will  be  an  important  area  of  study  of  the  DNA  HILAT 
satellite  mission  (Fremouw  et  al.,  1983).  The  current  convective  instability 
can  become  unstable  in  a  plasma  which  supports  a  density  gradient 
(perpendicular  to  the  ambient  magnetic  field)  and  a  current  which  flows 
parallel  to  the  ambient  3  field;  a  situation  which  can  occur  in  the  diffuse 
auroral  zone. 

The  Instability  was  Initially  studied  to  understand  plasma  disturbances 

in  laboratory  experiments  (Lehnert,  1958;  Hoh  and  Lehnert,  1960;  Kadomtsev  and 

Nedospasov,  1960)  but  has  been  more  vigorously  pursued  lately  in  regard  to 

ionospheric  disturbances.  The  linear  theory  of  the  mode  is  reasonably  well 

developed  In  the  short  wavelength  limit  (kL  »  1  where  k  is  the  wavenumber  and 

L  Is  the  scale  length  of  the  density  gradient).  Among  the  issues  considered 

thus  far  are  the  linear  properties  of  the  mode  in  the  low  altitude  auroral  F 

region  (Ossakow  and  Chaturvedi,  1979;  Vickrey  et  al. ,  1980),  the  high  altitude 

auroral  F  region  (Chaturvedi  and  Ossakow,  1981),  and  the  auroral  E  region 

(Chaturvedi  and  Ossakow,  1981);  as  well  as  studies  of  the  influence  of 

electromagnetic  effects  (Chaturvedi  and  Ossakow,  1981),  magnetic  shear  (Nuba 

and  Ossakow,  1980),  a  warm  electron  beam  (Chaturvedi  and  Ossakow,  1983), 
Manuscript  approved  November  9,  1983. 


.  velocity  shear  (Satyanarayana  and  Ossakow,  1983),  and  kinetic  effects  (s.g., 
finite  ion  Lannor  radius  effects,  wave-particle  resonances)  (Gary,  1983)  on 
the  instability.  A  nonlinear  theory  of  this  instability  has  also  been 
developed  (Chaturvedi  and  Ossakow,  1979)  and  numerical  simulations  of  the 
instability  have  been  performed  (Keskinen  et  al.,  1980). 

The  purpose  of  this  paper  is  to  present  analytical  and  numerical  results 
of  the  linear  properties  of  the  current  convective  instability  in  the  long 
wavelength  regime  (kL  <  1).  The  mathematical  analysis  is  similar  to  that  of 

the  long  wavelength  limit  of  the  E  x  3  instability  £Huba  and  Zalesak,  1983-). 

x  i nil  t  J  ^  y/ 

'ire'  derive^a  relatively  simple  dispersion  equation  of  the  mode,  and  presents 
simple  dispersion  relations  in  both  the  ion  coliisional  and  ion  inertial 
regimes.  These  analytical  results  are  compared  to  exact  numerical  results. 


The  organization  of  the  oaoer  is  as  follows. 


:r.e  next  section  we 


-derive — the  differential  equation  describing  the  mode.  la  SectriOC_li-± — we_ 

a-r-.L  )  3)  7  .  ’-ry  O' 

j?— present-  both  analytical  and  numerical  results.  Finally,  in  Section  IV  we 

tk>  r  o  /  J 

yi  -euoaarize -our- findings  and  discuss  applications  to  the  auroral  ionosphere. 


II.  DERIVATION  OF  MODE  EQUATION 

The  plasma  configuration  and  slab  geometry  used  in  the  analysis  are  shown 
in  Fig.  la.  The  ambient  magnetic  field  is  constant  and  in  the  z 
direction  { B  *  3^  e^)  ,  the  ambient  current  is  constant  and  in  the  z 

A 

direction  ( J  =  Jq  e^) ,  and  the  density  is  inhomogeneous  in  the  x  direction 

(n  *  ng(x))  .  A  weakly  coliisional  plasma  is  assumed  such  that  ve^/^e  <<  !■» 

v  /ft  «  1,  v.  /ft,  <<  l,and  v.  /ft  <<  1  (F  region  approximation)  where  ft  * 
en  6  16  x  xn  x  ct 

eB^/m^c  is  the  cyclotron  frequency  of  species  a,  refers  to  electron-ion 

collisions,  v  to  electron-neutral  collisions,  v,  to  ion-electron 

collisions,  and  v,  to  ion-neutral  collisions.  Furthermore,  we  asstime 

in 

that  v  /ft  <<  v.  /ft.  in  our  analysis.  Perturbed  quantities  are  assumed  to 
en  e  in  i 


vary  as  $p  ■  <5p(x)  exp  [  i( k^y  +  'a^z  -  oatl]  and  it  is  assumed  that  af?i.  «  1, 

kp ,  <<  1,  and  k  «  k  ,  where  p.  is  the  mean  ion  Larmor  radius.  That  is,  we 
i  z  /  i 

consider  low  frequency,  long  wavelength,  field-aligned  perturbations.  We  also 
assume  ,;Cz^\^pp  <<  1  where  k,jpp  is  the  electron  mean  free  path.  We  neglect 
temperature  effects.  Finally,  we  consider  only  electrostatic  oscillations  and 
assume  quasi-neutrality  (Qe  =  « 

The  fundamental  equations  used  in  the  analysis  are  continuity,  momentum 
transfer,  and  charge  conservation,  in  the  neutral  frame  of  reference: 
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where  a  denotes  species  (e:  electrons,  i:  ions)  and  other  variables  have 
their  usual  meaning.  Note  that  we  have  neglected  electron  inertia  effects  in 
Eq.  (2)  but  have  retained  ion  inertia  effects  in  Eq.  (3).  The  equilibrium 
drifts  are  given  by 
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where  u  *  u  -  kgVj  +  i»>  and  Che  prime  indicates  a  derivative  with  respect  to 
x.  We  have  neglected  collisional  effects  on  the  electron  motion  perpendicular 
to  B  and  on  the  ion  motion  parallel  to  B.  This  is  justified  since  we  have 
taken  ven/^e  <<  which  is  appropriate  for  auroral  ionospheric  conditions 
in  the  ?  region. 

We  now  substitute  Eqs.  (9)  and  (10)  into  Eq.  (4)  and  obtain 
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We  relate  5n  and  i  using  the  electron  continuity  equation  and  find  that  (from 
Eqs.  (1)  and  (9) 
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Finally,  substuciting  Eq.  (L2)  into  Eq.  (11)  we  arrive  ac  the  mode  equation 
for  the  current  convective  instability 
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III.  ANALYSIS  OF  MODE  EQUATION 

The  bulk  of  linear  analyses  of  the  current  convective  instability  have 

v  2 

made  use  of  the  local  approximation.  That  is,  it  is  assumed  that  k“L  » 
k“L  »  1  where  L  3  [3  b  n^/Sx)  is  the  density  gradient  scale  length  and  k.c 
*  3/3x.  With  this  assumption,  one  can  neglect  the  first  term  on  the  LHS  of 
Eq.  (13)  and  obtain  a  relatively  simple  dispersion  equation 


k2 


kV, 


k  a)  v  , 
y  ei 


.  i  k  ,  • .  ua 

,  2  ”  ca  '  k  ~  m  n0 

.  y  tu 


(14) 


which  has  been  thoroughly  analyzed  (Chaturvedi  and  Ossakow,  1981).  The  heart 
of  the  local  approximation  is  that  the  wavelengths  of  the  perturbations  are 
much  smaller  than  the  scale  length  of  the  density  gradient.  In  this  paper  we 
solve  Eq.  (13)  in  the  opposite  limit,  i.e.,  the  wavelengths  of  the 
perturbations  are  much  larger  than  the  scale  length  of  the  density  gradient 
(kL  «  1). 


We  consider  a  density  profile  with  a  single  discontinuity  at  x  =  0  (see 
Fig.  lb)  given  by 


nQ(x) 


x  >  0 


x  <  0 


For  x  *  0,  n^  *  0  and  Eq.  (13)  reduces  to 


where 


p  -  k“  T~p  -  0 
7  f 


«  k  a.fl  k  V 
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the  solution  to  Eq.  (16)  is  taken  to  be 


-k  Fx  k  Tx 
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Since  the  modes  are  assumed  to  be  bounded  as  x  +  +•  »  we  note  that 


-k  rx 

.  y 


x  >  0 


k  rx 

,  7 


x  <  0 


where  it  is  assumed  that  m  »  k^V^  (to  be  shown  a  posteriori). 

We  require  that  the  tangential  component  of  the  electric  field  be 
continuous  at  x  ■  0  (Hasegawa,  1971)  which  means  that  $  is  continuous  at 
x  *  0.  This  is  equivalent  to  requiring  that  the  interface  velocity  and  the 
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fluid  velocity  perpendicular  to  the  interface  be  equal  (Chandresekhar ,  19bl), 
i.a.,  5  is  continuous  at  the  discontinuity.  Thus,  6,  *  in  Eq.  (19)  so 
that 


-k  Tx 

v  y 

x  > 

0 

k  Tx 

*oe  7 

x  < 

0 

To  obtain  the  dispersion  equation  we  integrate  Eq.  (13)  across  the 
discontinuity  at  x  ■  0.  Thus,  we  have 


s  £  k  v  n . 

/  (nQr)dx  *  /  [n0kV*  -  Vo*-,dX 

— £  -£  Ui  3 


(21) 


Since  t>  is  continuous  across  the  boundary  at  x  =  0,  it  is  found  that  Eq.  (21) 
leads  to 


( V^I  "  L  V')2  *  K<ni  -  n2^o 

0) 


(22) 


where  (1,2)  indicate  the  region  x  >  0  (+s)  and  x  <  0  (-e),  respectively. 
Substituting  Eq.  (20)  into  Eq.  (22)  and  letting  £  +  0  ve  find  that 
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Equation  (23)  is  the  dispersion  equation  which  describes  the  long  wavelength 
modes  of  the  current  convective  instability.  From  Eq.  (23)  we  note 


(24) 
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and  u>  3  in  +  iv 


A.  Collisional  Limit 


We  consider  the  collisional  limit  v.  >>  id  so  that  cd 

in 


and  (24).  We  find  that 


n  L  n  -  n, 

_  >  v  -  —  — i _ ~ 

z  d  v .  r  n.  +  n_ 
ia  1  2 


ivia  in  Eqs.  (23) 
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Ln  ex 


where  the  subscript  ai  denotes  nonlocal.  Instability  can  occur  for  k2 7^ 

(nx  -  n2)/(nx  +  n,)  <  0. 

It  is  interesting  to  compare  the  growth  rate  of  the  instability  in  the 
long  wavelength  limit  (kL  «  1,  i.e.,  nonlocal)  to  that  of  the  short  wave¬ 
length  limit  (kL  »  1,  i.e.,  local).  Afsuming  v^n  >>  to,  we  find  from  Eq.  (14) 
that  the  short  wavelength  growth  rate  (local  growth  rate  y^)  is  given  by 


k  a.  v,  k  a,  a  . 

- E.  _ 1  f  1  +  _£_L_£_r1 

k  V,  L  (  ,  2  v.  v  , > 

y  in  k  in  ei 
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where  L  •  n  '/nn:  .  Defining  an  effective  wavenumber  k  »  [  k“  +  k“ 

^  0  0'  max  e  v  y  z 

1  /9 

lneni/vlnvei^  we  can  rewrite  Eq.  (25)  in  terms  of  Eq.  (27),  i.e., 
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is 


For  »  n->  we  have  Che  simple  relation  y^  =  kL  y  so  that  (1) 
proportional  to  the  ''wavenumber"  !<  and  (2)  <<  y ^  since  we  have  assumed  kL 

«  1. 

We  now  compare  the  analytical  expressions  derived  for  the  growth  rates 
(Eqs.  (25)  and  (27))  with  the  numerical  solution  of  Eq.  (13).  We  choose  a 
density  profile  given  by 


1  +  s  tanh(x/a) 
1  -  £ 


expression  (Eq.  (31))  in  both  the  short  wavelength  (kL  >>  1)  and  long 

wavelength  (k  <<  1)  limits.  Thus,  Eqs.  (25)  and  (27)  provide  good  estimates 
of  the  current  convective  instability  in  the  long  wavelength  and  short 
wavelength  limits,  respectively. 

B.  Inertial  Limit 

We  now  consider  the  ion  inertial  limit  given  bv  oi  >>  v .  so  that 

in 

m  m  in  Eqs.  (23)  and  (24).  The  dispersion  equation  is 
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We  obtain  a  simple  expression  for  the  growth  rate  by  assuming  that  « 

!’k"/k2'.  fij  /v  .  In  this  limit  Eq.  (32)  becomes 
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Again,  we  consider  the  normalization  m  -culL/Vjl  and  take  L/Vd  <  0.  We  find 
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In  this  limit  it  is  interesting  to  note  that  y  scales  as  (kyL)  and  is 
independent  of  kz/ky. 

We  contrast  Eq.  (34)  to  the  growth  rate  in  the  short  wavelength  limit  (kL 
»  1).  From  Eq.  (14)  we  find  that 
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Assuming  that  4  |  V^/Lft^  I  (  k^/k^j  /ftjjj  <<  1  we  find  that  Eq.  (36)  reduces  to 
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Following  Chaturvedi  and  Ossakow  (1981)  we  maximize  y ^  with  respect  to  k2/ky 
and  find  that  the  maximum  local  growth  rate  (y  )  in  the  inertial  limit  is 
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We  can  express  the  long  wavelength  growth  rate  in  terms  of  the  maximum  short 
wavelength  growth  rate  to  obtain 


r,  T  al  ~  n2 1  2/3 
Yni  aYlmL&y  n,  +  n»J 

j.  2 


where  a  is  a  numerical  factor  of  order  unity.  Thus,  for  n^  »  n2  we  find 
2/3 

that  y  „  »  (k  L)  y.  ,  in  contrast  to  the  collisional  limit  in 
tlx  7  xfli 

which  y^  scales  as  k^L. 

We  now  compare  the  analytical  expressions  for  the  growth  rate  with 

numerical  results  based  upon  Eq.  (13).  We  choose  the  same  density  profile  as 

-4 

in  the  collisional  case  (Eq.  (29)  with  e  »  0.8),  and  take  vg^/fte  *  10  , 


10"\  v.  /ft,  »  0,  and  |Lft,/V, |  »  104.  For  these  parameters,  the 
in  l  l  a 


assumption  that  led  to  Eq.  (36)  is  not  satisfied  and  the  analytical  expression 

a 

for  y ^  given  by  Eq.  (35)  must  be  used.  The  analytical  growth  rate  is  given  by 

0.62  k  L  »  1 

y 

Y  -  (39) 

0.86f k  Ll2/3  k  L  «  1 


The  results  of  the  comparison  are  shown  in  Fig.  3  where  we  plot  y  vs.  kvL. 

w 

A 

The  growth  rate  y  given  by  Eq.  (39)  is  plotted  as  shown  on  Fig.  (3)  and  the 
numerical  results  are  labelled  'exact.'  Again,  excellent  agreement  is  found 
between  the  analytical  and  numerical  results  in  both  the  short  wavelength  and 
long  wavelength  limits. 


IV.  DISCUSSION 

We  have  presented  an  analysis  of  the  current  convective  instability  in 

the  long  wavelength  limit  (kL  <<  1).  The  principal  result  of  the  paper  is  the 

derivation  of  a  relatively  simple  dispersion  equation  which  describes  the 

instability  in  the  long  wavelength  limit  (Eq.  (22)).  This  equation  is  solved 

analytically  in  two  limits:  collisional  (v,  >>  u)  and  inertial 

.  <<  tol.  We  find  that  in  the  collisional  limit  the  growth  rate  scales  as 

2/3 

k,  while  in  the  inertial  limit  it  scales  as  k  .  We  have  also  presented  a 
comparison  of  the  analytical  results  with  numerical  results  and  have  found 
very  good  agreement. 

We  now  discuss  the  application  of  these  results  to  the  auroral 
ionosphere.  Chaturvedi  and  Ossakow  (1981)  have  discussed  the  relevance  of  the 
short  wavelength  (or  local)  current  convective  instability  to  both  the  low 

altitude  (■*  400  km)  auroral  F  region  (collisional  limit)  and  to  the  high 

altitude  (■<•  1000  km)  auroral  F  region  (inertial  limit).  For  the  low  altitude 

-3  -1 

F  region  they  find  that  y  ^  ~  3  x  10  sec  ,  for  -  500 

2  -1  —2  —1 

m/sec,  v  .  '  5  x  10  sec  ,  v.  '  5  x  10  sec  ,  L  50  km  and  ma/m.>  ^ 
ex  in  e  x 

_  c  -.3 

3  x  10  ;  for  the  high  altitude  F  region  they  also  find  that  y  »  3  x  10 

-1  -1  -3  -1 

sec  but  for  30  sec  ,  <  10  sec  and  other  parameters  the  sane 
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as  the  low  altitude  F  region.  From  our  analysis  we  would  predict  that  the 
long  wavelength  modes  (kL  <  1)  would  have  much  longer  growth 

times  (1/y),  particularly  in  the  low  altitude  F  region. 

There  is  a  difficulty  in  applying  the  long  wavelength  results  to  the 

auroral  F  region  for  large  density  scale  lengths  (L  ■*  50  km).  The 

perpendicular  wavelengths  associated  with  the  long  wavelength  modes  are  such 
that  k,  L  <  1  which  leads  to  X,  >  2i  L  '  300  km.  However,  the  perpendicular 
spatial  scale  size  of  observed  auroral  F  region  density  enhancements,  which 
can  provide  the  zeroth  order  density  gradients  to  drive  the  instability,  can 
be  of  this  magnitude  (few  hundred  kms)  (Vickrey  et  al.,  1980;  Tsunoda  and 
Vickrey,  1983)  so  that  it  is  difficult  to  satisfy  the  condition  k^L  <  1.  On 
the  other  hand,  much  smaller  scale  sizes  of  auroral  structure  occur  during 
discrete  aurora.  Discrete  aurora  appear  to  have  two  distinct  scale  sizes. 

One  is  tens  of  kilometers  and  is  associated  with  inverted  V  precipitation. 
The  other  is  of  the  order  of  a  kilometer  and  is  associated  with  discrete 

auroral  arc  elements  (Davis,  1978).  The  auroral  arc  elements  appear  to  be 
imbedded  in  the  larger  inverted  V  structure.  Thus,  application  of  our  theory 
to  structure  in  discrete  auroral  arc  elements  leads  to  perpendicular 

wavelengths  A^  >  few  kms  which  can  satisfy  the  requirement  k, L  <  1. 

Finally,  several  aspects  of  the  present  theory  of  the  current  convective 

instability  deserve  comment.  First,  the  parallel  wavelengths  of  the 

irregularities  are  much  larger  than  the  perpendicular  wavelengths.  For 

-4  4 

typical  parameters  we  find  that  k^/k^  -  10  so  that  A^  -  10  A^.  The 
instability  in  the  long  wavelength  limit  can  then  produce  parallel  wave 

4 

structures  A^  >  10  km  which  is  larger  than  the  parallel  system  size.  Thus, 
it  is  important  to  develop  an  appropriate  theory  to  consider  the  finite 
parallel  extent  of  the  auroral  ionosphere  (e.g.,  Sperling  (1983)).  Second, 
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electromagnetic  effects  also  need  to  be  considered  in  a  more  detailed 
theory.  Chaturvedi  and  Ossakow  (1981)  investigated  these  effects  on  the  local 


theory  of  the  currenc  convective  instability  and  found  them  to  be 


negligible. 


However,  the  electromagnetic  corrections  are  proportional 


2  2  ? 

to  a)  / c  k"  and  may  become  large  in  the  long  wavelength  limit  since  kL  <  1 


We  are  presently  investigating  both  of  the  above  mentioned  effects. 


Acknowledgments 


We  thank  Mike  Keskinen  for  helpful  conversations  and  for  a  critical 
reading  of  the  manuscript.  This  research  has  been  supported  by  the  Defense 
Nuclear  Agency  and  the  Office  of  Naval  Research. 


Plasma  geometry  and  3lab  configuration  used  in  the  analysi 
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'EXACT 


Fig.  Z  Plot  of  the  growth  rate  [y  *  ylL/Vjl)  vs.  wavenumber  (k^L)  in  the 

collisional  limit  »  u>)  for  both  analytical  and  numerical 

-4  -4 

results.  The  parameters  used  are  v  , /ft  «  10  ,  v.  /ft.  -  10  ,  and 

ei  e  in  i  ’ 

labelled  accordingly.  The  numerical  results  are  based  upon  Eqs.  (13) 
and  (29)  with  e  *  0.8. 


r-- 0  62 


r  =  0.86(kyL); 


EXACT 


ig-  3  Plot  of  the  growth  rate  (y  *  ylL/V^I)  vs.  wavenumber  (k^L)  in  the 

inertial  limit  (v^n  <<  m)  for  both  analytical  and  numerical 

results.  The  parameters  used  are  v  . /Q  *  10  ,  k  /k  ■  10  ,  and 

ei  e  z  y  ’ 

*  0.  The  analytical  results  are  based  upon  Eq.  (39)  and 
labelled  accordingly.  The  numerical  results  are  based  upon  Eqs.  (13) 
and  (29)  with  e  ■  0.8. 
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DAHLGREN  LABORATORY 
DAHLGREN,  7A  22448 

OICY  ATTN  CODE  DF-L4  R.  3UTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
OICY  ATTN  COnE  465 
OICY  ATTN  CODE  461 
OICY  ATTN  CODE  402 
Old  ATTN  CODE  420 
OICY  ATTN  CODE  421 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

Old  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/ XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  30912 

Old  ATTN  SPDQQ 

Old  ATTN  X? 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 

Old  ATTN  OPR  HAROLD  GARDNER 
Old  ATTN  LXB  KENNETH  S.W.  CHAMPION 
OICY  ATTN  OPR  ALVA  T.  STAIR 

Old  ATTN  PHD  JURGEN  3UCHAU 

Old  ATTN  PHD  JOHN  ?.  MULLEN 

AF  WEAPONS  LABORATORY 
CRTLAND  AFT,  NM  37117 
Old  ATTN  SUL 

Old  ATTN  CA  ARTHUR  IT.  GUENTHER 
Old  ATTN  NTYCE  1LT.  G.  XRAJXI 

.AFT  AC 

PATRICK  .AFB,  FL  32925 
Old  ATTN  TF/MAJ  WILEY 
Old  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT- PATTERSON  AFB,  OH  45433 
Old  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  4  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
Old  .ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OP  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
Old  ATTN  J.  DBAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/ YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
Old  ATTN  YSEA 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OP  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

01CY  ATTN  DCKC  MAJ  JVC.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT -PATTERSON  AF3,  OH  45433 
01CY  ATTN  NICD  LIBRARY 

01CY  ATTN  ETDP  3.  BALLARD 


DEPARTMENT  OP  ENERCY 
AL3UQUERQUE  OPERATIONS  OFPICS 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG4G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

01CY  ATTN  DOC  CON  FOR  J.  3REEDL0  VE 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIPPISS  AF8 ,  NY  13441 

01CT  ATTN  DOC  LI3RARY/TSLD 
01CY  ATTN  OCSE  7.  COYNE 

SAMSO/SZ 

POST  OFFICE  30X  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  SZJ 

STRATEGIC  .AIR  COMMAND /IPFS 
OFTUTT  AFB,  MB  53113 

01CY  ATTN  ADWATE  MAC  3RUCS  3 AUER 
01CY  ATTN  MRT 

01CY  ATTN  DOK  CHIS?  SCIENTIST 

SAMSO/SK 
P.O.  30X  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CIA  VIM 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERCY 
LI3RARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOVITZ 


UNI7ERSITT  OP  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECa  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGEP. 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  37545 


OICY 

ATTN 

DOC 

CON 

FOR 

WOLCOTT 

OICY 

ATTN 

DOC 

COW 

FOR 

X, 

F.  TASCHEX 

OICY 

ATTN 

DOC 

CON 

FOR 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

r 

3  INN 

OICY 

ATTN 

DOC 

CON 

FOR 

?• 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTER  VE LT 

OICY 

ATTN 

D,  SAPPENFIELD 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  37115 

OICY  ATTN  DOC  CCN  FOR  W.  3R0WN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

QICT  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DI 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  3.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERCY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  TO  SONG 


OTHER  GO VEHEMENT 

DEPARTMENT  OP  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

OICT  (ALL  CORRES:  ATTN  SEC  OFFICER  FOR) 

INST ITU'  "R  TELECOM  SCIENCES 
NATIONAL  ^COMMUNICATIONS  4  INFO  ADMIN 
BOULDER,  CO  80303 

OICT  ATTN  A.  JEAN  (UNCLASS  ONLT) 

OICT  ATTN  W.  UTLAUT 
OICT  ATTN  D.  CR0MBI5 
OICT  ATTN  L.  3ERRT 

NATIONAL  OCEANIC  4  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICT  ATTN  R.  GRUBB 
OICT  ATTN  AERONOMT  LAB  C.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES ,  CA  90009 
OICT  ATTN  I.  GARFUNKEL 
01 CT  ATTN  T.  SALMI 
OICT  ATTN  V.  JOSEPHSGN 
OICT  ATTN  S.  3 OVER 
OICT  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01B03 

OICT  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
OICT  ATTN  L.  SLOAN 
OICT  ATTN  R.  THOMPSON 

3ERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  30X  983 
BERKELEY,  CA  94701 
OICT  ATTN  J. .  WORKMAN 
OICT  ATTN  C.  PRETTIE 
OICT  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  30X  3707 
SEATTLE,  WA  98124 
OICT  ATTN  C.  SISTER 
OICT  ATTN  D.  MURRAY 
OICT  ATTN  C.  HALL 
OICT  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICT  ATTN  D.3.  COX 
OICT  ATTN  J.P.  GILMORE 

C CMS AT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICT  ATTN  C.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NT  14850 

OICT  ATTN  D.T.  FARLEY ,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICT  ATTN  H.  LCC3T0N 
OICT  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

506  Wllshlra  31 vd. 

Santa  Monica,  Calif  90401 
OICT  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICT  ATTN  J.  ROBERTS 
OICT  ATTN  JAMES  MARSHALL 

CENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORCE  SPACE  CENTER 
GODDARD  BLVD  XINC  OF  PRUSSIA 
P.O.  30X  3555 
PHILADELPHIA,  PA  '.9101 

OICT  ATTN  M.H.  30RTNER  SPACE  SCI  LAB 

CENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NT  13201 
OICT  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC 
HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AR  99701 

(ALL  CLASS  AXTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  3R0WN  (UNCLASS  ONLY) 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  31730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  51820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CT  ATTN  K.  YEH 


JOHNS  HOPKINS  UNIVERSITY 
■APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL.  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  WARREN  S.  KNAPP 
01CT  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  30X  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 


INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  3EAURECARD  STREET 
ALEXANDRIA,  VA  22311 
Old  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
Old  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  3ENGST0N 


LOCKHEED  MISSILES  i  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

Old  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
01CT  ATTN  RICHARD  G.  JOHNSON  DEPT  52-12 
Old  ATTN  J.B.  CLAD  IS  DEPT  52-12 


INTL  TEL  &  TELECRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLET,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

l 101 l  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIECO,  CA  92138 
Old  ATTN  J.L.  SPERLING 


MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  30X  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
Old  ATTN  L.  LOUCHLIM 
Old  ATTN  D.  CLARK 
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MCDONNEL  DOUCLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  3 EACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
Old  ATTN  J.  HOULE 
Old  AITN  GEORGE  HR02 
Old  ATTN  3.  OLSON 
01CT  ATTN  R.W.  3ALPRIN 
Old  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATS  STREET 

SANTA  3AEL3ARA,  CA  93101 


Old 

ATTN 

P.  FISCHER 

Old 

ATTN 

W.?.  ORE TIER 

Old 

ATTN 

STEVEN  L.  GUTSCHZ 

Old 

ATTN 

R.  BOGUSCH 

Old 

ATTN 

R.  HENDRICK 

Old 

ATTN 

RALPH  XXL3 

Old 

ATTN 

DAVE  SOWLE 

Old 

ATTN 

F •  FA JEN 

01CY 

ATTN 

M.  SCHEI3E 

Old 

ATTN 

CONRAD  L.  LONCMIRE 

Old 

ATTN 

3.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
AL3U0UERQUE,  NEW  MEXICO  37106 
OICT  R.  5T2LLINCWERF 
Old  M.  ALMS 
Old  L.  VRICHT 

MITRE  CORPORATION,  THE 
?.0.  30X  208 
3EDF0RD,  MA  01730 

01CT  ATTN  JOHN  MORGANSTERN 
Old  ATTN  C.  HARDING 
Old  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  OOLLT  MADISON  3LVD 
MCLEAN,  VA  22101 
Old  ATTN  W.  HALL 
Old  ATTN  3.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  3LVD. 

LOS  ANCSLES,  CA  90025 

Old  ATTN  E.C.  FIELD,  JR. 


PENN  STL  VANIA  STATE  UNI  VERS  ITT 
IONOSPHERE  RESEARCH  LAB 
313  ELEdRICAL  ENGINEERING  EAST 
UNIVERSITT  PARK,  ?A  16302 
(NO  CLASS  TO  THIS  ADDRESS) 

Old  ATTN  IONOSPHERIC  RESEARCH  LA3 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

Old  ATTN  IRVING  L.  KCFSXT 

PHTSICAL  3TNAMICS,  INC. 

P.O.  BOX  3027 
3 ELLS VUE,  WA  98009 

Old  ATTN  S.J.  FREMOUW 

PHTSICAL  DYNAMICS ,  INC. 

P.O.  30X  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  S  D  ASSOCIATES 
P.O.  SOX  9695 
MARINA  DEL  RET,  CA  90291 
Old  ATTN  FORREST  GILMORE 
01CT  ATTN  WILLIAM  3.  WRICHT,  JR. 
01CT  ATTN  ROBERT  ?.  LSLS7TER 
01CT  ATTN  WILLIAM  J.  XAA2AS 
Old  ATTN  H.  ORY 
01CT  ATTN  C.  MACDONALD 
Old  ATTN  R.  TURCO 
01CT  ATTN  L.  D«RAND 
Old  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
Old  ATTN  CULLEN  CRAIN 
01CT  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  3ARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  VEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 
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SCIENCE  .APPLICATIONS ,  ISC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

Old  ATTN  LEWIS  M.  LINSON 
01CY  ATT"  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  ST RAKER 
01CT  ATTN  CURTIS  A.  SMITH 
Old  ATTN  JACK  MCDODGALL 

SCIENCE  APPLICATIONS,  INC 
17X0  GOODRIDGE  DR. 

MCLEAN,  7A  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


01CY 

ATTN 

DONALD  MEILSON 

01CY 

ATTN 

.ALAN  BURNS 

01CY 

ATTN 

G.  SMITH 

01CY 

ATTN 

R.  TSUNODA 

01CY 

ATTN 

DAVID  A.  JOHNSON 

01CY 

ATTN 

WALTER  G.  CHESNUT 

01CY 

ATTN 

CHARLES  L.  RINO 

01CY 

ATTN 

WALTER  JAYE 

01CY 

ATTN 

J.  VICKREY 

01CY 

ATTN 

RAY  L.  LEADA3RAND 

Old 

ATTN 

G.  CARPENTER 

Old 

ATTN 

G.  PRICE 

01CY 

ATTN 

R.  LIVINGSTON 

Old 

ATTN 

V.  GONZALES 

Old 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  .ATTN  W.P.  BOQUIST 


TRW  DEFENSE  S  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  3 EACH,  CA  90273 
01CY  ATTN  R.  K.  PLESUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 
01CY  ATTN  D/  STOCKWELL 
SNTF/1575 


7ISIDYNE 

SOUTH  3EDF0RD  STREET 
BURLINGTON,  MASS  01803 
01CT  ATTN  W.  SEIDY 
01CT  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 


TO YON  RESEARCH  CO. 

P.O.  Bo*  6890 
SANTA  3ARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 
01CY  ATTN  JOEL  GARBARINO 
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IONOSPHERIC  MODELING  distribution  LIST 
(UNCLASSIFIED  ONLY) 

PLEASE  DISTRIBUTE  ONE  COPT  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 
NOTED) 


NAVAL  RESEARCH  LA30RAT0RT 

NASA 

WASHINGTON,  D.C.  20375 

GODDARD  SPACE  FLIC3T  CENTER 

Dr.  P.  MANGE  -  CODE  4101 

GREENBELT,  MD  20771 

Dr.  ?.  COODMAN  -  CODE  4130 

DR.  K.  MAEDA 

A.F.  GEOPHYSICS  LABORATORY 

DR.  S.  CURTIS 

DR.  M.  DUBIM 

L.G.  HANSCCM  FIELD 

DR.  N.  MAYNARD  -  CODS  596 

BEDFORD,  MA  01730 

DR.  T.  E LEINS 

COMMANDER 

DR.  W.  S WIDER 

NAVAL  AIR  SYSTEMS  COMMAND 

MRS.  R.  SAGALYN 

DEPARTMENT  OF  THE  NAVY 

DR.  J.M.  FORBES 

WASHINGTON,  D.C.  20360 

DR.  T.J.  XENESKEA 

DR.  T.  CEUBA 

DR.  W.  3 CTREE 

OR.  H.  CARLSON 

COMMANDER 

DR-  J.  JASPERSE 

NAVAL  OCEAN  SYSTEMS  CENTER 

BOSTON  UNIVERSITY 

SAN  DISCO,  CA  92152 

MR.  R.  ROSE  -  CODS  5321 

DEPARTIENT  CF  ASTRONOMY 

BOSTON,  MA  02215 

NOAA 

DR.  J-  AARONS 

DIRECTOR  OF  SPACE  AND 

CORNELL  UNIVERSITY 

ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  30302 

ITHACA,  NY  14850 

DR.  A.  CLENN  JEAN 

DR.  W.E.  SWARTE 

DR.  C.W.  ADAMS 

DR.  D.  FARLEY 

DR.  D.N.  ANDERSON 

DR.  M.  KELLEY 

DR.  K.  DAVIES 

HARVARD  UNIVERSITY 

DR.  R.F.  DONNELLY 

HARVARD  SQUARE 

OFFICE  OF  NAVAL  RESEARCH 

CAMBRIDGE,  MA  02133 

800  NORTH  QUINCY  STREET 

DR.  M.B.  ttcELROY 

ARLINGTON,  VA  22217 

DR.  R.  LINDZSM 

DR.  C.  JOINER 

INSTITUTE  FOR  DEFENSE  ANALYSIS 

PENNSYLVANIA  STATE  UNIVERSITY 

400  ARMY/ NAVY  DRIVE 

UNIVERSITY  PARK,  PA  16802 

ARLINGTON,  VA  22202 

DR.  J.S.  NIS3ET 

DR.  E.  BAUER 

DR.  P.R.  ROHRBAUCH 

MASSACHUSETTS  INSTITUTE  OF 

DR.  L.A.  CARPENTER 

DR.  M.  LEE 

TECHNOLOGY 

DR.  R.  DIYANY 

PLASMA  FUSION  CENTER 

DR.  ?.  3ENNETT 

LIBRARY,  NU16-262 

DR.  F.  KLEVANS 

CAM3RIDCE,  MA  02139 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 

LA  JOLLA,  CA  92037 

DR.  D.A.  riAHLia 

DR.  E.  PRIEMAN 
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STANFORD  UNIVERSITY 
STANFORD ,  CA  94305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 

ballistic  research  laboratory 

A3EP.DSEN,  MD 
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